Stationary scalar configurations around extremal charged black holes 
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We consider the minimally coupled Klein-Gordon equation for a charged, massive scalar field in 
the non-extremal Reissner-Nordstrom background. Performing a frequency domain analysis, using 
a continued fraction method, we compute the frequencies cj for quasi-bound states. We observe 
that, as the extremal limit for both the background and the field is approached, the real part of 
the quasi-bound states frequencies lZ(u) tends to the mass of the field and the imaginary part 
T(ui) tends to zero, for any angular momentum quantum number I. The limiting frequencies in 
this double extremal limit are shown to correspond to a distribution of extremal scalar particles, at 
stationary positions, in no-force equilibrium configurations with the background. Thus, generically, 
these stationary scalar configurations are regular at the event horizon. If, on the other hand, the 
distribution contains scalar particles at the horizon, the configuration becomes irregular therein, in 
agreement with no hair theorems for the corresponding Einstein-Maxwell-scalar field system. 
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I. INTRODUCTION 

Scalar test fields in black hole (BH) geometries do 
not admit, generically, stationary configurations with an 
asymptotic decay and with real frequencies, i.e bound 
states. This follows from the physical requirement that 
only ingoing waves can exist at the horizon, therefore 
preventing a real equilibrium configuration between the 
field and the BH. Consequently, the configurations al- 
lowed in BH backgrounds are quasi-bound states, for 
which the frequencies are complex, with the imaginary 
part revealing a time dependence for the states, signal- 
ing either their absorption or, in the case of superradiant 
instabilities, their amplification by the BH 

The profile and some of the physical properties of 
quasi-bound states diverge at the horizon. This is inti- 
mately related to the inability that BH backgrounds have 
to accommodate, in a regular fashion, the scalar field, as 
an exact stationary solution, a property established by 
no-hair theorems 0, El] ■ But even with this caveat such 
quasi-bound states are informative. For instance in 
performing numerical simulations and starting with reg- 
ular initial data for a scalar field around a Schwarzschild 
BH, there were found damped oscillating solutions with 
frequency and decay rate described by the real and imagi- 
nary parts of quasi-bound state frequencies. These decay 
rates can be very small [f| and thus long lived scalar field 
configurations could exist around BHs, even though eter- 
nal and regular configurations are, in general, precluded 
by no-hair theorems. 

In this work we provide an example in which scalar field 
configurations around a BH can become stationary and 
an interpretation to justify why this is possible. We start 
by finding the quasi-bound states for a charged massive 



scalar field in a non-extremal Reissner-Nordstrom (RN) 
BH, and show that these states have a simple limiting 
behaviour as extremality for both the BH and the field is 
approached: the imaginary part of the frequency vanishes 
and the real part of the frequency tends to the test field 
mass (and charge). Then, the scalar field configurations 
obtained in this limit are understood as the electrostatic 
potential of a distribution of extremal scalar particles in 
equilibrium with the extremal BH, thus providing exam- 
ples of scalar fields around extremal charged BHs which 
are regular on the event horizon. Such configurations do 
not preserve spherical symmetry around the BH and that 
is the way they circumvents no-hair theorems. 

This paper is organized as follows. In Sec. |H]wc discuss 
solutions of the minimally coupled Klein-Gordon equa- 
tion on the RN background. Quasi-bound state solutions 
of this equation are considered in more detail in Sec. IIII1 
where some explicit frequencies are computed and ana- 
lyzed. In Sec. IIVI the extremal limit is discussed, by 
computing the states with frequencies equal to the limit- 
ing behaviour observed in Sec. IIII1 and an interpretation 
for these states is given. We draw some concluding re- 
marks and comment on the non-linear solution including 
the backreaction of the scalar field in Sec. [Vj 



II. BACKGROUND AND TEST FIELD 

We consider a massive, charged scalar field, $, with 
mass (i and charge q, obeying the wave equation 



D"D U 



$ = 



(1) 



where D u = D v — iqA v . This field is propagating in the 
background of a Reissner-Nordstrom BH with charge Q 
and mass M: 
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where f(r) = (r- r+)(r - r_)/r 2 , r± = M±y/AP - Q 2 
and A v dx v = —Q/rdt. Taking the standard ansatz for 
the scalar field which reflects the spherical symmetry and 
staticity of the background: 



* = 5>7* = £ e-^Y e m (e, <t>)R e (r) 
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where Y/™ are the spherical harmonics and w the complex 
frequency of a wave, ([T]) yields the radial equation for 
each mode: 
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where U = r 2 [(ojr - qQ) 2 - f{pL 2 r 2 + l(i + 1))] . Ob- 
serve that the azimuthal quantum number m is irrelevant 
due to spherical symmetry. The solution for the field $ 
is immediately obtained by solving the radial equation 
for each mode Q , due to the linearity of the wave equa- 
tion (JT|). In terms of Z(r) = rR(r), and dropping the 
subscript £ for notation simplicity, becomes 
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Vef f (r)] Z(r)=0,(5) 



where we have defined the effective potential 
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Alternatively, introducing the Regge- Wheeler radial co- 
ordinate r* , by dr* = dr/f(r), this wave equation is 
rewritten as 



dr* 2 



V, 



eff\ 



Z(r) = io 2 Z(r) , 



(7) 



where r — r(r*). The properties of the potential V e ff (r) 
have been discussed in the past, see eg. @. In particular 
one can show that the height of the centrifugal barrier 
increases with the charge of the field and that the con- 
stant value of the potential near the outer horizon also 
increases with the charge of the field but only up to some 
maximum; then it starts decreasing. The main feature of 
this potential, however, is that for a given combination of 
the parameters it exhibits a well, that can be considered 
as one of the key ingredients to have quasi-bound states. 

In order to solve the differential equation ([5]) we must 
provide a set of suitable boundary conditions at the hori- 
zon and at spacial infinity. To see the most relevant fea- 
ture of the near horizon behaviour we note that in this 
region equation ([7]) becomes to leading order: 



dr* 2 
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where </>+ = Q/r + is the electrostatic potential of the ex- 
ternal horizon. This equation is solved by a superposition 
of in and outgoing waves. Choosing the solution 



Z(r) 



(9) 



where w c = <?</>+, corresponding to an ingoing wave for 
q = 0, one observes the salient feature that it becomes 
an outgoing wave for lo < lu c (in this electromagnetic 
gauge). This is the condition for superradiance. 

Asymptotically, keeping the terms of order 1/r in equa- 
tion (jU) one gets 



R(r) 



where 
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(10) 



(11) 



Prom pop one observes a qualitatively distinct behaviour 
depending on the sign of the real part of x, H(x)- I n P ar_ 
ticular, for TZ(x) < we have quasi-bound states. These 
are characterized by a decaying behaviour at spatial in- 
finity. For 7Z(x) > we have scattering states. Hereafter 
we will be interested in quasi-bound states. 



III. SEMI-ANALYTIC GLOBAL SOLUTION: 
QUASI-BOUND STATES FREQUENCIES 

To find the solution of equation (|5|) in the region 
r > r + we will use a continued-fraction procedure de- 
veloped by Leaver to find the quasinormal modes for the 
Schwarzschild and Kerr BHs pJEl This amounts to take 
a power series ansatz with a pre-factor adapted to the 
boundary conditions observed in the previous section 
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where 
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Substituting (fT2"j) into (|5|) we obtain a three term recur- 
rence relationship for the a n of the form: 



a ai + P a a = , 



(14) 



1,2,3, 



with 



On = (1 - Q 2 )n 2 + c n - (1 - Q 2 ) + c , 
A* = -2(1 - Q 2 )n 2 + cm + c 2 , 
7„ = (1 - Q 2 )n 2 + c 3 n + c 4 . 

The constants Ci are lengthy expressions but otherwise 
straightforward to obtain. For a non charged massive 



1 I n [Hi Leaver considers a similar method for a Regge- Wheeler 
equation obtained from perturbing the Reissner-Nordstrom back- 
ground. 
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scalar field on a Schwarzschild background, these expres- 
sions reduce to the ones (in the same limit) given in Ref. 
Jj|. The procedure to find the frequencies of the quasi- 
bound states consists in setting the three term relation- 
ship (fT4)) in the form of a continued fraction algebraic 
equation and then solving it with a root finding proce- 
dure. 

Furthermore, we have checked the value of the frequen- 
cies obtained with the continued fraction method by nu- 
merically integrating the radial equation. We took as 
initial condition the behaviour of the function close to 
the horizon and integrated up to some large r compared 
with the horizon radius. 

The observed trend for the frequency of the quasi- 
bounded states is that, as the test field becomes extremal 
(p = |g|), the imaginary part of the frequency decreases. 
The imaginary part of the frequency is a measure of the 
rate at which the field falls into the horizon. For the 
charged massive scalar field, the closer the background 
black hole is to extremality the smaller this imaginary 
part is for the fundamental tone. The real part of the 
frequency, on the other hand, tends to increase and con- 
verge towards the field mass as the BH charge is increased 
to extremality This trend is shared by modes with differ- 
ent angular momentum quantum number. To make clear 
the behaviour of both the real and imaginary parts of the 
quasi-bound states frequencies as the double extremality 
is attained (|Q|, \q\ — > M,p), we plot them against the 
BH charge for q/p = 1 in Fig. [TJ Concerning the radial 
function as the double extremal limit is taken, the pro- 
file of the Ro(r) mode is qualitatively different from the 
profiles for the others modes. When £ = (s-wave), the 
radial function is localized in the region between the ex- 
ternal horizon and the maximum of the potential barrier. 
As \Q\/M — > 1, Ro(r) becomes narrower in such a way 
that its maximum tends to the horizon. For £ ^ 0, the 
functions Ri(r) tend to spread out as the double extremal 
limit is approached - Fig [5J 



IV. EXTREMAL BLACK HOLE 

The results of the previous section describe the lim- 
iting behaviour of the quasi-bound states as the double 
extremal limit is attained. We shall now consider ex- 
actly this limit by focusing on the extremal Reissner- 
Nordstrom BH, |Q| = M (r± = M), with an extremal 
test field (p = \q\). Then, using a new radial coordinate 
p = r — M, for states with u> = p, = \q\ and qQ > 0, the 
radial wave equation ^ reduces to 



d 2 R e 2 dR e l{l + 1) 
dp 2 p dp p 2 



R f = 



(15) 



This equation is the radial part of the Laplace equation 
on Euclidean 3-space E 3 whose solution is Rg = Agp 1 + 
Bi/ p +1 . The spacial part of the scalar field <f> is therefore 
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FIG. 1. (Top panel) The imaginary part of the frequency 
always tends to zero, even if the rates of decay do depend on 
the value of /i. (Bottom panel) The real part of the frequency 
tends to /i as \Q\/M —¥ 1. In both plots we keep fi/\q\ = 1, 
£ = 1 and M = 1. 



a linear combination of harmonic functions 

B f 



-ifit 



A e p 



(16) 



Each of these partial waves, with appropriate Ag, Bp, de- 
scribes the double extremal limit of a quasi-bound state. 

To interpret the meaning of the modes (fTB) and un- 
derstand their appearance, it is useful to rewrite the ex- 
tremal RN background using the coordinate p; this cor- 
responds to isotropic coordinates. Then the fields take 
the form 



ds 2 



-H~ z dt 



H 2 5 l3 dx l dx 3 



A = H- X dt , (17) 



where furthermore p = ^/ 1 b~ijX % x3 = |x| and H is a har- 
monic function on Euclidean 3-space with a simple pole 
localised at the origin: H = 1 + M/|x|. In these coor- 
dinates x = is the location of the extremal RN BH 
horizon. 

Taking the scalar field in the form $(t, x l ) = 
e~ ttJ,t H(x l ), the wave equation ([T]) with p = \q\ in the 
background (fl"7|) , yields the harmonic equation A^H = 
0. One solution is the harmonic function with a single 
pole at x': H = /i/|x — x'|. This describes the electric 
potential of one particle located at this pole. Expressed 
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FIG. 2. The real part of the radial function for some values of 
the BH charge. We have also plotted the effective potential for 
Q — 0.7 and q = p, = 0.4. (Top panel) For I — the maximum 
of the radial function is approaching the outer horizon, as 
extremality of the background is approached. (Bottom panel) 
For 1=1 the function tends to spread along the potential well. 



in terms of spherical coordinates on E 3 , (p,9,4>), chosen 
such that x' lies at coordinates (p\ 9' ^ 0, (/>'), then 



H = 



1 l,m 



B,(p',9',^') 



where 



4-7T/J, 

2£+ 1 



(18) 



(19) 



For the latter, a pure s-wave corresponds to a source 
localised at the horizon, whereas in the former, the peak 
of the radial function tends to the origin, as displayed in 
the top panel of Fig. [2] A related observation is that the 
harmonic function H is regular at the horizon except if 
the source is localised there. 

A similar interpretation for the modes in (|16[) carries 
through if H represents multiple scalar sources instead of 
a single one. We then have a superposition of harmonic 
functions with localised poles at fixed positions x^., cor- 
responding to spherical coordinates (p' k , 9' k , tfi' k ), 



i 



(20) 



which may again be rewritten, in spherical coordinates, 
as the multipolar expansion in the right hand side of (|18l) , 
with Bi(p' ,9' ,(f)') replaced by J2k Be(p'k,9'k, ^4)> corre- 
sponding to replacing one particle by many particles. 

The existence of stationary scalar states and their in- 
terpretation can, furthermore, be generalized to a back- 
ground with multiple extremal BHs instead of a single 
one. This is achieved replacing H by a superposition 
of harmonic functions with localised poles at different 
points, Xj. H = 1 + Y^i Mi/\x — Xj|, since the scalar field 
equation still reduces to a harmonic equation on E 3 . Such 
solution of the Einstein-Maxwell system is the well known 
Majumdar-Papapetrou multi BH solution [13, HH, corre- 
sponding to a collection of BHs with mass and charge 
Mi = Qi, placed at arbitrary positions Xj [HI, held 
in equilibrium by a balance between gravitational and 
electrostatic forces. A multiple scalar particle configura- 
tion will be regular on each horizon of the Majumdar- 
Papapetrou background as long as Xj ^= x' fc , for all i, k. 

The scalar particles are in equilibrium with the BHs 
due to a 'no-force' condition, a balance between grav- 
itational and electromagnetic forces, as can be easily 
checked by studying the orbits generated by the La- 
grangian C = [i\J 

— 9avi a i v + pA a x a in the background 
([XT]) , where 'dot' denotes derivative with respect to 
proper time. This Lagrangian is adequate to describe the 
interaction of the scalar particles with the background 
(|17[) because, at linear level, there is no interaction me- 
diated by the scalar field; only gravitational and electro- 
magnetic interactions occur. The gravitational energy 
added to the system in equilibrium - the multi BH so- 
lution - by the massive scalar field is balanced by the 
electromagnetic energy carried by the field. 



So H is indeed the spatial part of (fT6[) . with Ai = 0. 
This fact shows that the radial profiles (|16p correspond 
to partial waves for an extremal point-like scalar source 
of mass and charge p displaced from the BH horizon, i.e. 
at p 7^ 0. If, one the other hand, the particle is at p = 0, 
it follows from ([T8 )) - (fT9]) that only the s-wave appears. 

There is an agreement between the behaviour of the 
s-wave seen in Sec. IHII as extremality is approached, and 
the one exhibited for the extremal case in this section. 



V. CONCLUSIONS 

A scalar field on a BH geometry does not, generically, 
admit stationary configurations. In this note we have 
showed that an extremal scalar field in the background 
of a charged, extremal BH geometry does admit such 
configurations and we have provided a physical interpre- 
tation for them. 
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Our first observation was that the frequencies of quasi- 
bound states of a massive, charged, minimally coupled 
scalar field in the RN background have a well defined 
behaviour when a double extremal limit, for both the 
test field and the background is taken: the imaginary 
part vanishes and the real part becomes equal to the 
field mass. Then we showed that in such double extremal 
limit, configurations with a real frequency equal to the 
particle's mass exist, corresponding to a distribution of 
extremal scalar particles, placed at arbitrary locations in 
the exterior of the extremal (multi-)BH solution. If none 
of these particles sits at the BH horizon, the configuration 
is regular therein. One may argue, however, that the field 
is irregular at the location of the sources. But this is the 
traditional problem in classical field theory associated to 
point-like sources. 

The stationary scalar field states we have exhibited arc 
due to no-force configurations between scalar sources and 
extremal BHs, at linear level: the gravitational attraction 
is being balanced by the electromagnetic repulsion. At 



non-linear level, however, the scalar field will back react 
on the geometry and, since it is charged, it will source 
the Maxwell field wherever the scalar field is non-trivial 
and not just at the location of the sources, in contrast 
to the typical multi-centre solutions found in Supergrav- 
ity/String theory (see, e.g. 13]). It would will be inter- 
esting, but also challenging, to study the configurations 
we have analyzed herein at non-linear level, as solutions 
of the corresponding Einstcin-Maxwcll-scalar field the- 
ory. 
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